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Abstract 

In this paper, we investigate the power efficient resource allocation algorithm design for secure 
multiuser wireless communication systems employing a full-duplex (FD) base station (BS) for serving 
multiple half-duplex (HD) downlink (DL) and uplink (UL) users simultaneously. We propose a multi¬ 
objective optimization framework to study two conflicting yet desirable design objectives, i.e., total 
DL transmit power minimization and total UL transmit power minimization. To this end, the weighed 
Tchebycheff method is adopted to formulate the resource allocation algorithm design as a multi¬ 
objective optimization problem (MOOP). The considered MOOP takes into account the quality-of- 
service (QoS) requirements of all legitimate users for guaranteeing secure DL and UL transmission 
in the presence of potential eavesdroppers. Thereby, secure UL transmission is enabled by the FD BS 
and would not be possible with an HD BS. The imperfectness of the channel state information of 
the eavesdropping channels and the inter-user interference channels is incorporated for robust resource 
allocation algorithm design. Although the considered MOOP is non-convex, we solve it optimally by 
semidefinite programming (SDP) relaxation. Simulation results not only unveil the trade-off between 
the total DL transmit power and the total UL transmit power, but also confirm the robustness of the 
proposed algorithm against potential eavesdroppers. 
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I. Introduction 

The exponential growth in high data rate communieation has triggered a tremendous demand 
for radio resourees sueh as bandwidth and energy. An important technique for reducing the 
energy and bandwidth consumption of wireless systems while satisfying quality-of-service (QoS) 
requirements is multiple-input multiple-output (MIMO), as it offers extra degrees of freedom for 
efficient resource allocation. However, the MIMO gain may be difficult to achieve in practice due 
to the high computational complexity of MIMO receivers. As an alternative, multiuser MIMO 
(MU-MIMO) has been proposed as an effective technique for realizing the MIMO performance 
gain. In particular, in MU-MIMO systems, a transmitter equipped with multiple antennas (e.g. a 
base station (BS)) serves multiple single-antenna users which shifts the computational complexity 
from the receivers to the transmitter O. Yet, the spectral resource is still underutilized even if 
MU-MIMO is employed as long as the BS operates in the traditional half-duplex (HD) mode, 
where uplink (UL) and downlink (DL) communication are separated orthogonally in either time 
or frequency which leads to a significant loss in spectral efficiency. 

Full-duplex (FD) wireless communication has recently received significant attention from both 
academia and industry due to its potential to double the spectral efficiency of the existing 
wireless communication systems Q-lIHl. In contrast to conventional HD transmission, FD enables 
simultaneous DL and UL transmission at the same frequency. However, in practice, a major 
challenge in FD communication is the self-interference (SI) caused by the signal leakage from 
the DL transmission to the UL signal reception. Although Q, dH reported that SI can be 
partially cancelled through analog circuits and digital signal processing, the residual SI still 
severely degrades the performance of FD systems if it is not properly controlled. Besides, co¬ 
channel interference (CCI) caused by the UL transmission impairs the DL transmission. Thus, 
different resource allocation designs for FD systems were proposed and studied to overcome 
these challenges. For example, the authors of |l5l investigated the end-to-end outage probability 
of MIMO FD single-user relaying systems. In d^, a resource allocation algorithm was proposed 
for the maximization of the end-to-end system data rate of multi-carrier MIMO FD relaying 
systems. In d3, massive MIMO was applied in FD relaying systems to facilitate SI suppression 
and to improve spectral efficiency. Simultaneous DL and UL transmission via an FD BS in small 
cells was studied in dll, where a suboptimal DL beamformer was designed to improve the system 
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throughput. 

On the other hand, seeurity is a erueial issue for wireless eommunieation due to the broadcast 
nature of the wireless medium. Traditionally, secure communication is achieved by cryptographic 
encryption performed at the application layer and is based on the assumption of limited compu¬ 
tational capabilities of the eavesdroppers. However, new computing technologies (e.g. quantum 
computers) may make this assumption invalid which results in a potential vulnerability of tradi¬ 
tional approaches to secure communication. The pioneering work in dH proposed an alternative 
approach for providing perfectly secure communication by utilizing the nature of the channel in 
the physical layer. Specifically, |l9l revealed that secure communication can be achieved whenever 
the information receiver enjoys better channel conditions than the eavesdropper. Inspired by this 
finding, multiple-antenna transmission has been proposed to ensure communication security ifTOll - 
ini, since multiple antennas provide spatial degrees of freedom which can be utilized to degrade 
the eavesdropper channel. In particular, transmitting artificial noise (AN) is an effective means to 
deliberately impair the information reception at the eavesdroppers ifTOll . In ifTTI . a power allocation 
algorithm was designed for maximizing the secrecy outage capacity via AN generation. The 
authors of ifT^ proposed a transmit beamforming approach for secrecy provisioning by generating 
spatially selective AN. In lfT3ll . joint transmit signal and AN covariance matrix optimization was 
studied for secrecy rate maximization. However, all of the above works focused on HD systems 
and the obtained results may not be applicable to FD communication systems. In fact, for FD 
communication systems, both UL and DL users are exposed to the risk of eavesdropping because 
of the simultaneous UL and DL transmission. Therefore, it is necessary to ensure communication 
security for DL and UL concurrently. Although guaranteeing DL security with a multiple-antenna 
HD BS has been exhaustively studied in the literature [[TT]| - [[T3ll . UL communication cannot be 
secured with an HD BS which can perform either transmission or reception in each time instant 
but not both, and thus, cannot jam the eavesdroppers in the UL. On the other hand, the single¬ 
antenna UL users lack the required spatial degrees of freedom to ensure communication secrecy. 
Multiple-antenna FD BSs are a promising solution to this problem due to their inherent capability 
of performing simultaneous transmission and reception. 

The notion of secure communication in FD systems has received some attention recently. In 
HU, joint information beamforming and jamming beamforming for an FD BS was proposed 
to guarantee DL and UL security. Yet, l[T4l assumed that there is no CCI between DL and UL 
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users and that the SI at the FD BS ean be eaneelled perfeetly, whieh may be too optimistie 
assumptions for praetieal FD systems. Besides, lfT4ll also assumed that the eavesdropper ehannel 
state information (CSI) was perfeetly known at the FD BS whieh whieh is highly idealistie. 
In faet, some idle users (e.g. roaming users) in the system may misbehave and eavesdrop the 
information signal of the legitimate users. Perfeet CSI of these potential eavesdroppers may not 
be available due to their diseontinuous interaetion with the BS. The authors of ifTSll proposed 
an optimal power alloeation algorithm to guarantee seeure eommunieation for an FD reeeiver 
employing only statistieal CSI of the eavesdropper ehannel. In [|T^ . robust beamforming for 
the ease of imperfeet CSI was studied for two-way FD eommunieation systems. However, the 
seeure eommunieation approaehes proposed in ifTSll and ifT^ eannot be applied in FD MU-MIMO 
wireless eommunieation systems direetly due to the eonsiderable differenees in the eonsidered 
system models. Furthermore, to the best of the authors’ knowledge, power effieient seeure FD 
eommunieation has not been studied in the literature yet. Speeifieally, total DL and total UL trans¬ 
mit power minimization are eonflieting design objeetives in seeure FD eommunieation networks. 
In our previous work lH], we investigated a power effieient resouree alloeation algorithm for FD 
systems under a multi-objeetive optimization framework whieh unveiled a trade-off between total 
DL and total UL power eonsumption. However, studying this trade-off is still an open problem 
for seeure FD systems. Besides, in |[T1, perfeet CSI knowledge was assumed. 

In this paper, we address the above issues. To this end, the resouree alloeation algorithm 
design for seeure FD eommunieation networks is formulated as a multi-objeetive optimization 
problem (MOOP). The proposed MOOP formulation jointly minimizes the total DL transmit 
power and the total UL transmit power for seeure MU-MIMO wireless eommunieation systems 
employing an FD BS for guaranteeing both UL and DL seeurity. Besides, our problem formu¬ 
lation takes into aeeount the imperfeetness of the CSI of the links between the FD BS and the 
potential eavesdroppers, the links between the UL users and potential eavesdroppers, and the 
CCI links. Although the eonsidered MOOP is non-eonvex, we solve it optimally by semidefinite 
programming (SDP) relaxation leading to a set of Pareto optimal resource allocation policies. 
Our simulation results not only unveil the trade-off between the total DL transmit power and 
the total UL transmit power, but also confirm the robustness of the proposed algorithm against 
imperfect CSI. 
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II. System Model 

In this section, we present the considered MU-MIMO FD wireless communication system 
model. 

A. Notation 

We use boldface capital and lower case letters to denote matrices and vectors, respectively. 
A^, Tr(A), Rank(A), and det(A) denotes the Hermitian transpose, trace, rank, and determinant 
of matrix A, respectively; A”^ and A^ represent the inverse and Moore-Penrose pseudoinverse 
of matrix A, respectively; A ^ 0, A 0, and A ^ 0 indicate that A is a positive semidefinite, 
a positive definite, and a negative semidefinite matrix, respectively; Ia? is the N x N identity 
matrix; denotes the set of all A x M matrices with complex entries; EI^ denotes the 

set of all A X A Hermitian matrices; l-j, H-H, and H-HiT’ denote the absolute value of a complex 
scalar, the Euclidean vector norm, and the Frobenius matrix norm, respectively; £{■} denotes 
statistical expectation; diag(a;i, • • • ,xk) denotes a diagonal matrix with the diagonal elements 
given by {xi, • • • , xk} and diag(X) returns a diagonal matrix having the main diagonal elements 
of X on its main diagonal. extracts the real part of a complex-valued input; [x]"*" stands for 
max{0,x}; the circularly symmetric complex Gaussian distribution with mean // and variance 
cr^ is denoted by CAA(yU., cr^); and ~ stands for “distributed as”. 

B. Multiuser System Model 

We consider a multiuser communication system. The system consists of an FD BS, K legiti¬ 
mate DL users, J legitimate UL users, and M roaming users, cf. Figure [IJ The FD BS is equipped 
with At > 1 antennas for facilitating simultaneous DL transmission and UL reception in the 
same frequency bancj^. The K + J legitimate users are single-antenna HD mobile communication 
devices to ensure low hardware complexity. The number of antennas at the FD BS is assumed 
to be larger than the number of UL users to facilitate reliable UL signal detection, i.e.. At > J. 
Besides, the DL and the UL users are scheduled for simultaneous UL and DL transmission. 
Unlike the local legitimate signal-antenna users, the M roaming users are travelling wireless 
devices from other communication systems and are equipped with Ar > 1 antennas. The 

'We note that circulator based FD radio prototypes, which can transmit and receive signals simultaneously on the same 
antennas, have been demonstrated El- 
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Fig. 1. A multiuser communication system with a full-duplex (FD) radio base station (BS), K = 1 half-duplex (FID) DL users, 
J = 1 HD UL users, and M = 1 HD roaming user (potential eavesdropper). 

multiple-antenna roaming users are searehing for aeeess to loeal wireless servicej^ However, 
it is possible that the roaming users deliberately intereept the information signal intended for 
the legitimate users if they are in the same service area. As a result, the roaming users are 
potential eavesdroppers which have to be taken into account for resource allocation algorithm 
design to guarantee communication security. In this paper, we refer to roaming users as potential 
eavesdroppers and we assume A^t > A^r for studying resource allocation algorithm design. 

C. Channel Model 

We focus on a frequency flat fading channel. In each scheduling time slot, the FD BS transmits 
K independent signal streams simultaneously at the same frequency to the K DL users. In 
particular, the information signal to DL user k E {1,..., K} can be expressed as 

Xfc = (1) 

where E C and E are the information bearing signal for DL user k and 

the corresponding beamforming vector, respectively. Without loss of generality, we assume 

However, the signal intended for the desired user may be eavesdropped by the roaming users. 
Hence, in order to ensure secure communication, the FD BS also transmits AN to interfere the 
reception of the roaming users (potential eavesdroppers). Therefore, the transmit signal vector, 

^In order to receive the wireless services provided by the local FD BS, the roaming users have to transmit pilot signals to 
facilitate system clock synchronization and channel estimation. 
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X G comprising K information streams and AN, is given by 

K 

x = ^Xfc + z, (2) 

k=l 


where z G represents the AN veetor generated by the FD BS to degrade the ehannel 

quality of potential eavesdroppers. In partieular, z is modeled as a eomplex Gaussian random 
veetor with z ~ CAA(0, Z), where Z G ZAO, denotes the eovarianee matrix of the AN. 
Therefore, the reeeived signals at DL user /c G {1, ..., AT}, the FD BS, and potential eavesdropper 
m G {1,..., M} are given by 

K j 




,UL 


= hf Xfc 



h^z 


+ 




■UL 

■j 




i=i 


y = 


artifieial noise 

multiuser interferenee eo-ehannel interferenee 

J K 

+ Hsi^Xfc + and 

V ^ artifieial noise 

self-interferenee 


( 4 ) 


K 


J 


= (5) 

artifieial noise 

respeetively. The DL ehannel between the FD BS and user k is denoted by G and 

fj k G C represents the ehannel between UL user j and DL user k. Variables 8{\d^^\^} = 1, 
and Pj are the data and transmit power sent from UL user j to the FD BS, respeetively. Veetor 
gj G C^txi denotes the ehannel between UL user j and the FD BS. Matrix Hsi G C^tx^t 
denotes the SI ehannel of the FD BS. Matrix G denotes the ehannel between the 

FD BS and potential eavesdropper m. Veetor ^ denotes the ehannel between UL 

user j and potential eavesdropper m. Variables h^, /j fc, gj, Hsi, L^, and „j eapture the joint 
effeet of path loss and small seale fading, ~ CAf(0, ~ CAA(0, and 

nj^ ~ CAA(0, (tI^Iatj^) represent the additive white Gaussian noise (AWGN) at the FD BS, DL 
user k, and potential eavesdropper m, respeetively. In ([3]), the term J2j=i V^fj,kd^^ denotes the 
aggregated CCI eaused by the UL users to DL user k. In ([4]), the term Hsi Ylk=i represents 


the SI. 
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III. Resource Allocation Problem Formulation 

In this section, we first define the adopted performance metrics for the considered multiuser 
communication system. Then, we discuss the assumptions regarding the CSI knowledge for 
resource allocation. Finally, we formulate the resource allocation problems for DL and UL 
transmit power minimization, respectively. For the sake of notational simplicity, we define the 
following variables: A: e {1,..., K}, and Gj = gjgf, j G {1,..., J}. 


A. Achievable Rate and Secrecy Rate 


Assuming perfect CSI at the receiver, the achievable rate (bit/s/Hz) of DL user k is given by 


= iog.,(i + r 


DL\ 
k ) 


with = 


h?wJ 


E |hf w.P + EEI/.aP + Tr(H,Z) + 
i=i 


( 6 ) 


r^k 




where F^^ is the receive signal-to-interference-plus-noise ratio (SINK) at DL user k. Besides, 
the achievable rate of UL user j is given by 


R 


UL 


= iog2(i + r 


UL^ 

j > 


with 


Pjigfv 


pUL ^ 

^ j 


X:Pn|g^v,f + Tr(pV,diag (HsiZHg+ f Hgiw.wfHg 

n^j ^ ^ k=l 


( 7 ) 

( 8 ) 


+0-1 


ULl 


where F^^ is the receive SINK of UL user j at the FD BS. The variable G is 

the receive beamforming vector for decoding the information received from UL user j and we 
define Vj = G In this paper, zero-forcing receive beamforming (ZF-BF) 

is adopted. In this context, we note that ZF-BF closely approaches the performance of optimal 
minimum mean square error beamforming (MMSE-BF) when the noise term is not dominating^ 
ifTTll or the number of antennas is sufficiently large [|71- Besides, ZF-BF facilitates the design of 
a computational efficient resource allocation algorithm. Hence, the receive beamformer for UL 
user j is chosen as = (u^Q^)^, where = [O,..., 0,1,0,..., O], Q'l’ = (Q^Q)"^Q^, and 

(i-i) (J-j) 

Q = [gi,..., gj]. The term Tr (^pVj diag (^HsiZH|^j -h Ef=i Hgi^Vfc^vf Hg j j in dH]) models 
the impact of the imperfectness of the SI cancellation [fT^ Eq. (4)] due to the limited receiver 
dynamic range and 0 < p <C 1 is a constant modelling the noisiness of the SI cancellation at the 


^We note that the noise power at the BS is not expected to be the dominating factor for the system performance since BSs 
are usually equipped with a high quality low-noise amplifier (LNA). 
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FD BS. In particular, lfT9]| has shown that this model accurately captures the combined effects 
of additive automatic gain control noise, non-linearities in the analog-to-digital converters and 
the gain control, and oscillator phase noise which are present in practical hardware. 

As outlined before, for guaranteeing communication security, roaming users are treated as 
potential eavesdroppers who eavesdrop the information signals desired for all DL and UL 
users. Thereby, we design the resource allocation algorithm under a worst-case assumption for 
guaranteeing communication secrecy. In particular, we assume that a potential eavesdropper can 
cancel all multiuser interference before decoding the information of a desired user. Thus, under 
this assumption, the channel capacity between the FD BS and potential eavesdropper m for 
eavesdropping desired DL user k and the channel capacity between the UL user j and potential 
eavesdropper m for overhearing UL user j can be written as 

= log 2 det(lM, + X-‘L"wiwfL„) and (9) 

= log2det(lA,, + F,.X-‘e,,„e«„), (10) 


respectively, where X„i = L^ZL^ + denotes the interference-plus-noise covariance 


matrix for potential eavesdropper m. We emphasize that, unlike an HD BS, the FD BS can 
guarantee both DL security and UL security simultaneously via AN transmission. The achievable 
secrecy rates between the FD BS and DL user k and UL user j are given by 


j^DL—Sec 
^UL-Sec 


— max I C, 


DL-E 

k.m 


and 


R 


UL 


max I C, 


UL-El 

J 


( 11 ) 

( 12 ) 


respectively. 

B. Channel State Information 


In this paper, we focus on slowly time-varying channels. At the beginning of each time slot, 
the FD BS obtains the CSI of all channels to facilitate global resource allocation. In practice, 
the UL users perform handshaking with the FD BS which facilitates UL channel estimation at 
the FD BS. Since the channels may change slowly in time, the UL users embed pilot signals 
periodically in the data packets. Hence, the FD BS is able to frequently update and refine the 
CSI estimate of the UL users. Furthermore, for the acquisition of the CSI of the DL users at the 
FD BS, handshaking is also performed between the FD BS and the DL users at the beginning 
of each scheduling slot which allows the FD BS to obtain the statuses and QoS requirements of 
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the DL users. Then during transmission, the DL users are required to send acknowledge (ACK) 
packets to inform the FD BS of successful reception of the data packets. Hence, the FD BS 
can regularly update the CSI estimates of the DL transmission links. Therefore, perfect CSI 
for the UL and DL transmission links, i.e., gj, Vj G J}, and h^, V/c G K}, is 

assumed over the transmission period. On the other hand, for the CCI channels, the DL users 
can receive the pilot signal^ of the UL users and feed back the CCI channel estimates to the 
FD BS only at the beginning of each scheduling time slot. Hence, the FD BS can update the 
CSI of the CCI channels only at the beginning of every scheduling time slot. As a result, the 
CSI of the CCI channels at the FD BS is imperfect. The roaming users (potential eavesdroppers) 
also perform handshaking at the beginning of a scheduling slot which facilitates the estimation 
of the corresponding channels at the FD BS. However, the FD BS cannot update the channel 
information of the potential eavesdroppers during transmission since they are silent in the current 
time slot. Therefore, only imperfect CSI of the channels between the FD BS and the potential 
eavesdroppers is available. Furthermore, for the UL user-to-potential eavesdropper channels, 
although there is no direct interaction between the potential eavesdroppers and the UL users, the 
UL users can obtain the CSI by measuring the potential eavesdroppers’ pilot signals when the 
potential eavesdroppers perform handshaking with the FD BS. Then, the UL users can feed back 
the CSI of these channels to the FD BS. However, since the potential eavesdroppers only perform 
handshaking at the beginning of each slot, the FD BS can update the information of the UL user- 
to-potential eavesdropper channels only once in every scheduling time slot. Consequently, only 
imperfect CSI of the UL user-to-potential eavesdropper channels can be obtained at the FD BS. 
To capture the impact of imperfect CSI, we model the CSI uncertainty based on a deterministic 
model [l20]| - ll22]| . In particular, the CSI of the link between UL user j G {1,..., J} and DL user 
A; G {1,..., K}, i.e., fj^k, the CSI of the link between the FD BS and potential eavesdropper 
m G {1,..., M}, i.e., L^, and the CSI of the link between UL user j G {1,..., J} and potential 


''We assume that the DL users, UL users, and roaming users utilize orthogonal sequences as pilot signals, which allows the 
FD BS to distinguish the pilot signals of different users. 
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eavesdropper me {1,..., M}, i.e., are modeled as 

fj,k = fj,k + ^j,k — £ C : \Afj^k\ < (13) 

Lm = Lm + ALm, = |Lm ^ C^TxAfa . ||/\L^||^ < |, and (14) 

^j,m ^j,rn ^ ^ • ||Aej j^|| ^ ^ > (13) 

respeetively, where fc, L^, and are the CSI estimates available at DL user k, the FD 
BS, and UL user j, at the beginning of a scheduling slot, respectively; Afj^k, AL^, and Aoj^m 
denote the unknown channel uncertainties due to the time varying nature of the channel. The 
continuous sets and contain all possible channel uncertainties with bounded 

magnitude ej^k, SDLm^ and £ul,,„., respectively. In practice, the values of Sj^k, £dl,„, and £ul,,„. 
depend on the coherence time of the associated channels and the transmission duration of the 
scheduling slot. 


C. Optimization Problem Formulation 


We first study the problem formulation for two desirable system design objectives of the 
considered secure FD communication system. Then, we investigate the two system design ob¬ 
jectives jointly under a multi-objective optimization framework. The first considered objective is 
the minimization of the total DL transmit power at the FD BS and is given by 


Problem 1 (Total DL Transmit Power Minimization): 

K 

'2 '.Tr(Z) 


minimize A ||’Wfc| 


s.t. Cl: min 




J 




E | hfw ,|2 + ^ F ^.|/,.,|2 + Tr ( H , Z ) 


a- 


r^k 


i=i 


life 


JT,, 12 


C2: 




> V) 

— ^ req.- 5 ^ J't 




EPn|g^v,f+ Tr pV^diag HsiZHg+ EHsiw^wfHi^ 


K 


k=l 


req. 


+^ULl|Vi|P 


C3: max logs det(IjvR + L^) < i?toi, V/c, m, 

ALiTTiGifiDLyTi ’ 

C4: niax logs det {In^ + ^, Vj, m, 

AL/TTiEifiDLyn ’ 

Aej^eriuLj_„, 

C5: Pj > 0, Vi, C6: Z ^ 0. 


(16) 


The system design objective in (fT^ is to minimize the total DL transmit power which is 
comprised of the DL signal power and the AN power. Constants > 0 and > 0 
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in constraints Cl and C2 in (fT^ are the minimum required SINK for DL users k E {1,... ^K} 
and UL users j E respeetively. In partieular, in eonstraint Cl, the minimum required 

SINK for DL user k is satisfied for a given CSI uneertainty set ^ for the CCI ehannels. 

, in C3 and C4, respeetively, are pre-defined system parameters representing 
the maximum tolerable data rate a^iotential eavesdropper m for deeoding the information of DL 
user k and UL user j, respeetiveljo In faet, DL and UL seeurity is guaranteed by eonstraints C3 
and C4 for given CSI uneertainty sets and OuLj rn ■ partieular, if the above optimization 

problem is feasible, the proposed problem formulation guarantees that the seereey rate for DL 
user k is bounded below by > log 2 (l + r°q^) — and the seereey rate 

for UL user j is bounded below by > iog,(i + rs-1 — }. Constraint C5 

777 , ^Jy'^ 

is the non-negative power eonstraint for UL user j. Constraint C6 and Z E are imposed 
sinee eovarianee matrix Z has to be a Hermitian positive semidefinite matrix. We note that the 
objeetive of Problem 1 is to minimize the total DL transmit power under eonstraints C1-C6 
without regard for the eonsumed UL transmit powers. 

The seeond system design objeetive is the minimization of total UL transmit power and ean 
be mathematieally formulated as 

Problem 2 (Total UL Transmit Power Minimization): 

J 

minimize y Pj 

ZGH^T,wfe,P,- 

j=l 

s.t. Cl - C6. (17) 


Problem 2 targets only the minimization of the total UL transmit power under eonstraints C1-C6 
without taking into aeeount the total eonsumed DL transmit power. 

The objeetives of Problems 1 and 2 are desirable for the system operator and the users, 
respeetively. However, in seeure FD wireless eommunieation systems, these objeetives eonfiiet 
with eaeh other. On the one hand, the DL information and AN transmission eause signifieant SI 
whieh impairs the UL signal reeeption. Henee, the UL users have to transmit with a higher power 
to eompensate this interferenee to satisfy the minimum required reeeive SINK of the UL users at 


^If the eavesdroppers do not emit pilot signals, the estimation errors for the eavesdropper channels, ALrn and Aej^rn, are 
random and follow certain distributions. In this case, the proposed resource allocation algorithm design can still be used but 
constraints C3 and C4 have to be converted to probabilistic constraints which specify the maximum tolerable secrecy outage 
probability [HI Eq. (30), (31)]. 
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the FD BS. On the other hand, a high UL transmit power results in a strong CCI for DL signal 
reeeption and a higher risk of information leakage to the potential eavesdroppers. Henee, the FD 
BS has to transmit both the DL information and the AN with higher power to ensure the QoS 
requirements of the DL users and the seeurity requirements of the DL and UL users. However, 
this in turn eauses high SI and gives rise to an esealating inerease in transmit power for both 
UL and DL transmission. To overeome this problem, we resort to multi-objective optimization 
[l2^ . If2^ . In the literature, multi-objective optimization is often adopted to study the trade-off 
between conflicting system design objectives via the concept of Pareto optimality ll2^ . Il24l . To 
facilitate our presentation, we denote the objective function of Problem i as Qj(wfc, Ti^Pj). The 
Pareto optimality of a resource allocation policy is defined in the following: 

Definition / l23l/ .- A resource allocation policy, {wfc,Z,Pj}, is Pareto optimal if and only if 
there does not exist any {w^, Z, Pj} with Qi(wfc, Z, Pfi) < Qfiwk, Z, Pj), Vi G {1, 2}. 

In other words, a resource allocation policy is Pareto optimal if there is no other policy that 
improves at least one of the objectives without detriment to the other objective. In order to capture 
the complete Pareto optimal set, we formulate a third optimization problem to investigate the 
trade-off between Problem 1 and Problem 2 by using the weighted Tchebycheff method ll2^ . 
The third problem formulation is given as 

Problem 3 (Multi-Objective Optimization): 



minimize max 

ZeH^T,wj,,Pj *=1>2 


s.t. Cl - C6 


(18) 


where Qi(wfc, Z, Pfi = -f Tr(Z) and Q 2 {^k, Z, Pj) = Pj- Q* is the optimal 


objective value of the z-th problem and is treated as a constant for Problem 3. Variable Aj > 0, 
^ • Aj = 1, specifies the priority of the Lth objective compared to the other objectives and 
reflects the preference of the system operator. By varying Aj, we can obtain a complete Pareto 
optimal set which corresponds to a set of resource allocation policies. Thus, the operator can 
select a proper resource allocation policy from the set of available policies. Compared to other 
formulation methods for handing MOOPs in the literature (e.g. the weighted product method 
and the exponentially weighted criterion (231), the weighted Tchebycheff method can achieve 
the complete Pareto optimal set with a lower computational complexity, despite the non-convexity 
(if any) of the considered problem. It is noted that Problem 3 is equivalent to Problem i when 
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Aj = 1 and Xj = 0, Vi 7 ^ j. Here, we mean by equivalenee that both problems have the same 
optimal solution. 


IV. Solution of the Optimization Problem 


Problems 1, 2, and 3 are non-convex problems due to the non-convex constraints C1-C4. 
Besides, constraints Cl, C3, and C4 involve infinitely many inequality constraints due to the 
continuity of the corresponding CSI uncertainty sets. To solve these problems efficiently, we 
first transform Cl, C3, and C4 into equivalent linear matrix inequality (LMI) constraints. Then, 
Problems 1,2, and 3 are solved by semidefinite programming (SDP) relaxation. 

To facilitate the SDP relaxation, we define W^. = and rewrite Problems 1-3 in the 

following equivalent forms: 

Equivalent Problem 1: 

K 

minimize TrlW^) + Tr(Z) 


s.t. Cl: min 


Tr(HfcW, 


g ^ + Tr(HfcZ) 

r^k j=l 


> 

— reqfc 






nfc 


C2: 


P,Tr(G,V, 


■>r 


UL 


EPn 


Sn^jl 


K 


+ Tr pV,- diag HsiZHg + ^ HsiW^H 


k=l 


— req,- 


+^ulI 


Vi, 


C3: 

C4: 


m^ log2det(l7VR + XjL^WfcL^) < VA;,m, 

fiDLm, ’ 

niax logs det(IiVR + ^, Wj, m, 

AL/TTiGiiDLyTi ’ 


C5: Pj > 0, Vi, C 6 : Z P 0, C7: P 0,VA:, C 8 : Rank(Wfc) < 1,VA:, 


(19) 


where W^. P 0, G and Rank(Wfc) < 1 in (fT9l) are imposed to guarantee that 

Wfc = holds after optimization. 

Equivalent Problem 2: ^ 

minimize y Pj 

Wfc,zeH^T,p,. ^ 

j=i 


s.t. C1-C8. 

Equivalent Problem 3: 


( 20 ) 


minimize r 

Wk,zm^T,Pj,T 

s.t. C1-C8, 

C9: A,(Q,-g*)<RV^e{l,2}, 


( 21 ) 
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where r is an auxiliary optimization variable and (| 2 TI) is the epigraph representation of (fTSl) . 

Sinee Problem 3 is a generalization of Problems 1 and 2, we foeus on solving Problem 3. 
Now, we introduee a Lemma whieh allows us to transform eonstraint Cl into an LMI. 

Lemma 1 (S-Procedure l^25\l }: Let a function fm{x),Tn E {l,2},x G be defined as 

/m(x) = x^A^x + 23fJ{b^x} + Cm, (22) 


where A^, G G and Cm G Then, the implication /i(x) < 0 ^ / 2 (x) <0 

holds if and only if there exists a variable <5 > 0 such that 


Ai 

bi 


A 2 

b2 

1 

cr 

Cl 


bf 

C2 


^ 0 , 


(23) 


provided that there exists a point x such that /^(x) < 0. 

To facilitate the presentation, we first define ffc = , fj^k] h = [fi,k, • • •, fj,k] Af^ = 

[A/i fc,..., A/jfc]^, and P = diag (Pi,... , Pj), where f^, ffc, and Affc denote the collections of 
the CCI channels, CCI channel estimates, and CCI estimation errors at DL user k, respectively. 
Hence, the collection of CCI channels at DL user k can be modeled as 


ffc = ffc 


Affc, G {1,..., K}, Dfc A <1 f, e : II Affcll < £fc k 


(24) 


where el = Ylij=i ^\k- As a result, by applying (l24l) . Cl can be equivalently expressed as 
Cl: 0>AffPAffc + 23f?{ffPAffc} + ffPffc + 5^Tr(HfcWk-^^^^^^ + Tr(HfcZ) + a^ 


r^k 


pDL 

reqj, 


2 


By exploiting Lemma 1, we obtain the following implications: 


Affc^Affc 




Cl 


holds if and only if there exists a variable ()fc > 0 such that 

U: R^^(Wfc,Z,P„4) 

4lj-P -ffc 

-ff -5fc4-<-ffc^Pffc-Tr(HfcZ 


B 


K 

Ew. 

r^k 


Wfc 

pDL 


Bh, P0,VA;,j, (25) 


holds, where Bhj, = [O hfc]. 

Next, for handling the non-convex constraints C3 and C4, we establish the following propo¬ 
sition for facilitating the constraint transformation for the considered optimization problem. 
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Proposition 1: For > 0 and > 0, we have the following implieations for 

eonstraints C3 and C4 of equivalent Problems 1-3, respeetively: 

C3 ^ C3: ^ VL^ G VAi,m, (26) 

C4 <(=» C4: Ve^-^ G f2uL,-„,, VL^ G Vj, m, (27) 

where = 2 — 1 and = 2 — 1. We note that C3 and C3 are equivalent 

respeetively if Rank(Wfc) < 1. Besides, C4 and C4 are always equivalent. 

Proof: Please refer to Appendix A. ■ 

Although C3 and C4 are eonvex LMI eonstraints whieh are less diffieult to handle eompared 
to C3 and C4, they still involve an infinite number of inequality eonstraints. To eireumvent this 
diffieulty, we introduee the following Lemma to further simplify C3 and C4. 

Lemma 2 (Generalized S-Procedure /|2^ ).~ Let /(X) = X^AX + X^B + B^X + C, and 
D ^ 0. For some f > 0, /(X) A 0, VX G |x| Tr(DXX^) < l|, is equivalent to 

A 0 . (28) 


and express eonstraint C3 as 


c 

B^ 


I 0 



— t 


B 

A 


- 1 

Q 

1 

o 


As a result, we first substitute L^ = L^ + AL^ into 


m 

0 :< AL"({gLz_W0AL„ + AL«(4Oi.Z-Wt)L. 

+ L2(4°i.Z - Wj)AL„ + L"(4°^Z - Wi)L„ + VI;, m, (29) 

for AL„ e |AL„|Tr(£,(? AL.„AL") < 1 i. Then, by applying Lemma 2, eonstraint C3 is 


equivalently represented as 

C3: ^ (Wfc, Z, tk^m) 


dllHz 










-2 


A^t 


-BL„,WfcBL„, P 0,VA;,m, (30) 


for tk^m > 0, yk,m, where Bl„^ = [Lm Iatt]- On the other hand, for eonstraint C4, two 
estimation error variables are involved, namely ^ and AL^, and Lemma 2 eannot be 
direetly applied. Henee, we introduee a slaek matrix variable Mj m G to handle the 

eoupled estimation error variables in eonstraint C4. In partieular, eonstraint C4 ean be equivalently 
represented by 

” ( 31 ) 


C4a. PjGj fYiGj ,^ A ^ ^ 7 )^) 

C4b: M,-™ A (^UL _ ^ 


(32) 
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Proposition 2: Constraint C4 holds if there exists a Hermitian matrix ^ G j G 

{1,..., J}, m G {1,..., M} whieh meets eonstraints C4a and C4b. 

Proof: Please refer to Appendix B. ■ 

Then, we apply Lemma 2 to eonstraints C4a and C4b in a similar manner as for handling C3 
and obtain the following equivalent LMIs for C4a and C4b, respeetively: 


C4a. Pjj 0!j,m) 

for aj^rn > 0, Vj, m, and 




Pj^j,m 


-P 


Pj T 


-2 


P0,yj,m, (33) 


C4b: 

l^ZL™ + _ M,,^ 


A0,Vj,m, (34) 


for > 0, Vj, m. 

The remaining non-eonvex eonstraint in (1211) is the rank-one eonstraint C 8 . Solving sueh a 
rank-eonstrained problem is known to be NP-hard [l27ll . Henee, to obtain a traetable solution, 
we relax eonstraint C 8 : Rank(Wfc) < 1 by removing it from the problem formulation, sueh that 
the eonsidered problem beeomes a eonvex SDP and is given by 

minimize r 

Wfe ,ZgH^t eH^R, 


s.t. C2,C5,C6,C7,C9, 


CIO. ^ 0, VA^, TTZ, 


Cl: R^^(Wfe,Z,P,,4) ^0,VA:, 

C4a. (]V[j „j, Tj, ttj A 0,Vj, 777 , 


C3: R^^_^(Wfc,Z, 4 ,^) A 0,Wk,m, 

C4b: Mj^rn, Pj,m) h 0,Vj,m. (35) 


The relaxed eonvex problem in (1351) ean be solved effieiently by standard eonvex program solvers 
sueh as CVX [|^ . Besides, if the solution obtained for a relaxed SDP problem is a rank-one 
matrix, i.e., Rank(Wfc) = 1 for 7 ^ 0, \/k, then it is also the optimal solution of the original 
problem. Next, we verify the tightness of the adopted SDP relaxation in the following theorem. 

Theorem 1: Assuming the eonsidered problem is feasible, for > 0, we ean always obtain 
or eonstruet a rank-one optimal matrix whieh is an optimal solution for (l35l) . 

Proof: Please refer to Appendix C. ■ 

By Theorem 1, the optimal beamforming veetor ean be reeovered from by performing 
eigenvalue deeomposition of and seleetion of the prineiple eigenveetor as w^. 






TABLE I 
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System parameters. 


Carrier center frequency and system bandwidth 

1.9 GHz and 200 kHz 

Path loss exponent and SI cancellation constant, p 

3.6 and -80 dB (3) 

DL user noise power and UL BS noise power, and a^L 

-100 dBm and -110 dBm^ 

Potential eavesdropper noise power, BS antenna gain 

— 100 dBm and 10 dBi 

Maximum tolerable data rate at potential eavesdroppers for DL users, fltoii, 

1 bit/s/Hz 

Maximum tolerable data rate at potential eavesdroppers for UL users, RITaj ^ 

1 bit/s/Hz 


V. Results 

In this section, we investigate the performance of the proposed multi-objective optimization 
based resource allocation scheme through simulations. The most important simulation parameters 
are specified in Table HI There are K = 3 DL users, J = 7 UL users, and M = 2 potential 
eavesdroppers in a cell. The users and potential eavesdroppers are randomly and uniformly 
distributed between the reference distance of 30 meters and the maximum service distance of 
600 meters. The FD BS is located at the center of the cell and equipped with iVr antennas. Each 
potential eavesdropper is equipped with TYr = 2 antennas. The small scale fading of the DL 
channels, UL channels, CCI channels, and eavesdropping channels are modeled as independent 
and identically distributed Rayleigh fading. The multipath fading coefficients of the SI channel 
are generated as independent and identically distributed Rician random variables with Rician 

factor 5 dB. To facilitate the presentation, we define the maximum normalized estimation error 

£2 

of the CCI channels, DL eavesdropping channels, and UL eavesdropping channels as = 

2 

||i °^||2 = 11 ^^, respectively. Lurthermore, we assume that different 

channels have the same maximum normalized estimation error, i.e., fTj f, = = ^Lt- 

In addition, we assume that all DL users and all UL users require the same minimum SINRs, 
respectively, i.e., T^^^ = and TU^^, = 

A. Transmit Power Trade-ojf Region 

In Ligure [2l we study the trade-off between the DL and the UL total transmit powers for a 
maximum normalized channel estimation error of = 5% and different numbers of antennas 
at the LD BS. The trade-off region is obtained by solving (1351) for different values of 0 < A* < 

®The DL user and potential eavesdropper noise powers are higher than the UL noise power at the BS since the BS is usually 
equipped with a high quality LNA 1291 . 
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Fig. 2. Average system objective trade-off region achieved by the proposed resource allocation scheme. The double-sided 
arrows indicate the power saving due to additional antennas. 

1, Vi G {1,2}, i.e., the A* are varied uniformly using a step size of 0.01 subjeet to = 1- We 
assume a minimum required DL SINK of = 10 dB and a minimum required UL SINK of 
Treq = 5 dB. It ean be observed from Figure [2] that the total UL transmit power is a monotonically 
deereasing funetion with respeet to the total DL transmit power. In other words, minimizing the 
total UL power eonsumption leads to a higher power eonsumption in the DL and viee versa. This 
result confirms that the minimization of the total UL transmit power and the total DL transmit 
power are conflicting design objectives. For the case of iVr = 12, 5 dB in UL transmit power can 
be saved by increasing the total DL transmit power by 6 dB. In addition. Figure 2 also indicates 
that a significant amount of transmit power can be saved in the FD system by increasing the 
number of BS antennas. This is due to the fact that the extra degrees freedom offered by the 
additional antennas facilitate a more power efficient resource allocation. However, due to channel 
hardening, there is a diminishing return in the power saving as the number of antennas at the 
FD BS increases ifTTll . 

For comparison, we consider a baseline scheme. For the baseline scheme, we adopt ZF-BF 
as DL transmission scheme such that the multiuser interference is avoided at the DL legitimate 
users. In particular, the direction of beamformer for legitimate user k is fixed and lays in the 
null space of the other DL legitimate users’ channels. Then, we jointly optimize Z, Pj, and the 
power allocated to under the MOOP formulation subject to the same constraints as in (l35l) via 
SDP relaxation. Figure |2] depicts the trade-off region for the baseline resource allocation scheme 
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for comparison. As can be seen, the trade-off regions achieved by the baseline scheme are above 
the curves for the proposed optimal scheme. This indicates that the proposed resource allocation 
scheme is more power efficient than the baseline scheme for both DL and UL transmission. 
Indeed, the proposed resource allocation scheme fully exploits the available degrees of freedom 
to perform globally optimal resource allocation. On the contrary, for the baseline scheme, the 
transmitter is incapable of fully utilizing the available degrees of freedom since the direction of 
the transmit beamformer is fixed. Specifically, the fixed beamformer can cause severe 
SI and increase the risk of eavesdropping which results in a high power consumption for UL 
transmission and the AN. Besides, the trade-off region for the baseline scheme is strictly smaller 
than that of the proposed optimal scheme. For instance, when At = 12, the baseline scheme 
can save only 1 dB of UL transmit power by increasing the total DL transmit power by 2 dB, 
due to the limited flexibility of the baseline scheme in handling the interference. 

We note that we also considered two other baseline schemes. For the first scheme, we adopted 
an isotropic radiation pattern for Z and optimized W^, and Pj. For the second scheme, we 
considered a wireless communication system with an HD BS. However, both of these schemes 
could not satisfy the QoS requirements in constraints Cl and C2 for the adopted simulation 
setting. Therefore, performance results for these two schemes are not shown in the paper. 

B. Average Total Transmit Power versus Minimum Required SINK 

In Figure [3l we investigate the average power consumption of the DL and UL users versus 
the minimum required DL SINK, Fj^^, for a minimum required UL SINK of = 5 dB and 
a maximum normalized channel estimation error of = 5%. The FD BS is equipped with 
At = 10 antennas. We select the resource allocation policy with Ai = 0.1 and A 2 = 0.9 which 
indicates that the system operator gives priority to total UL transmit power minimization. It 
can be observed that both the DL and the UL power consumption increase with F^q. However, 
the DL power consumption grows more rapidly than the UL power consumption. The reason 
behind this is threefold. First, a higher DL transmit power is required to fulfill the more stringent 
DL QoS requirements. Second, a higher AN transmit power is required to neutralize the higher 
potential for information leakage. In particular, with the growth of the DL transmit power, the 
channel capacity between the FD BS and the potential eavesdroppers increases. Hence, the FD 
BS also has to allocate more power to the AN to prevent interception by potential eavesdroppers 
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Fig. 4. Outage probability versus the minimum required 
DL SINR (dB), for different resource allocation 

schemes. 

which results in a higher DL power eonsumption. Third, the SI beeomes more severe for higher 
DL transmit powers. Specifically, since the FD BS tries to keep the total UL transmit power 
low to avoid strong CCI for inereasing required SINR, more degrees of freedom at the FD BS 
have to be utilized to suppress the SI to improve the receive SINR of the UL users. As a result, 
fewer degrees of freedom are available for reducing the DL transmit power consumption causing 
even higher DL transmit powers. On the other hand, the proposed seheme provides substantial 
power savings compared to the baseline seheme in both DL and UL for the entire range of F^g. 
Also, the baseline scheme eannot satisfy the QoS requirements when F^g is larger than 12 dB. 
In Figure HI we study the outage probability for the proposed seheme and the baseline seheme 
versus the minimum required DL SINR for a maximum normalized ehannel estimation error of 
^ resouree alloeation policy with Ai = 0.1 and A 2 = 0.9, and different numbers of 
antennas. An outage event oeeurs whenever the problem in (1351) is infeasible. As ean be observed, 
a large outage probability reduetion ean be aehieved with the proposed scheme compared to the 
baseline seheme. For the ease of F^g = 12 dB and At = 10, the outage probability for the 
proposed seheme is only 0.5%, whereas the outage probability for the baseline seheme is 99.3%. 
These results eonfirm that the proposed scheme is more robust and reliable in the presenee of 
imperfeet CSI eompared to the baseline seheme. 

In Figure [5l we study the power eonsumption of the DL and UL users versus the minimum 
required UL SINR, rj:^q, for a minimum required DL SINR of Fj^g = 10 dB and a maximum 




Fig. 3. Average power consumption (dBm) versus the 
minimum required DL SINR (dB), F^q, for different 
resource allocation schemes. 
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Fig. 5. Average power consumption (dBm) versus the minimum required UL SINK (dB), Fj^q, for different resource allocation 
schemes. 

normalized channel estimation error of = 5%. The FD BS is equipped with A^t = 10 
antennas. The other system parameters are identical to those adopted in Figure [H As can be 
observed, the total transmit power of DL and UL increase as the minimum required UL SINK 
increases. In fact, for the UL users, they need to transmit with higher power to fulfill the more 
stringent UL QoS requirements. Furthermore, increasing the total UL transmit power causes 
more CCI to the DL users. Hence, the FD BS has to allocate more power to the DL information 
signals to satisfy the DL QoS requirements. Moreover, the higher UL and DL information signal 
powers increase the susceptibility to eavesdropping for both UL and DL users. Hence, the FD 
BS also has to increase the AN power to guarantee DL and UL transmission security. On the 
other hand. Figure [5] clearly shows that the proposed scheme provides significant power savings 
compared to the baseline scheme in all considered scenarios. In fact, the baseline scheme cannot 
find a feasible solution for Fj^g > 8 dB. 

C. Average User Secrecy Rate versus Minimum Required SINR 

Figure [6] depicts the average user secrecy rate of the DL and UL users versus the minimum 
required DL SINR, F^g, for a maximum normalized channel estimation error of = 5% 
and = 5 dB. The FD BS is equipped with At = 10 antennas. We select the resource 
allocation policy with Ai = 0.1 and A 2 = 0.9. The average user secrecy rates for DL and UL 
users are calculated by averaging the total DL and the total UL secrecy rates, i.e., - 

J j^h — Sec 

and J -, respectively. It can been seen that the average DL user secrecy rate increases 
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[^Proposed scheme, DL, = 5% 

^Proposed scheme, UL, = 5% 

- B - Baseline, DL, = 5% 

- © - Baseline, UL, = 5% 

Minimum required DL user secrecy rate 

- V - Minimum required UL user secrecy rate 



6 8 10 12 14 16 18 2 

Minimum required DL SINR (dB) 

Fig. 6. Average user secrecy rate (bits/s/Hz) versus the 
minimum required DL SINR (dB), F^q. 


—3|l— Proposed scheme, DL, = 5% 


^ Proposed scheme, UL, = 5% 


- B - Baseline, DL, = 5% 


- 0 - Baseline, UL, = 5% 


Minimum required DL user secrecy rate 


-^-Minimum required UL user secrecy rate 


. i - 8 . il 




6 8 10 
Minimum required UL SINR (dB) 

Fig. 7. Average user secrecy rate (bits/s/Hz) versus the 
minimum required UL SINR (dB), Fjr^q. 


with r°q since the channel eapaeity between DL user k and potential eavesdropper m is limited 
to = 1 bit/s/Hz. Besides, the average UL user secrecy rate depends only weakly on 

due to the proposed robust optimization. In addition, we eompare the average DL and UL user 
seereey rates of the proposed seheme with the minimum required DL and UL user seereey rates, 
i.e., log 2 (l + and log 2 (l + respectively. As can be 

seen, although the CSI of the DL eavesdropping ehannels, UL eavesdropping ehannels, and CCI 
channels is imperfeet, the average user seereey rate achieved by the proposed seheme fulfills 
the minimum required user seereey rate in both DL and UL which confirms that the security 
of both links can be guaranteed simultaneously. This is due to the robustness of the proposed 
optimization algorithm design. On the other hand, the baseline seheme aehieves a higher seereey 
rate than the proposed seheme for r°q ranges from 2 dB to 12 dB. However, to accomplish this, 
the baseline scheme requires exceedingly large transmit powers at both the FD BS and the UL 
users eompared to the proposed scheme, cf. Figure [3l Besides, the superior performance of the 
baseline scheme in terms of seereey rate also eomes at the expense of a extremely high outage 
probability. In partieular, the baseline seheme is always infeasible when F^^ is larger than 12 
dB, ef. Figure IH 

Figure |7] illustrates the average user secrecy rate versus the minimum required UL SINR, Fj^q, 
for a maximum normalized ehannel estimation error of Kggt = 5% and a minimum required DL 
SINR of r^q = 10 dB. The system setting is the same as for Figure 0 It ean be seen that 
the average UL user seereey rate increases with F^q since the ehannel eapaeity between UL 
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Fig. 8. Average power consumption (dBm) versus the maximum normalized channel estimation error, 

user j and potential eavesdropper m is limited to = 1 bit/s/Hz. Besides, the average 

DL user seereey rate is not sensitive to the inerease of Similar to Figure l6l the proposed 
seheme achieves a higher average user secrecy rate than the minimum required user secrecy rate 
despite the imperfect CSI which confirms the robustness of the proposed scheme for guaranteeing 
communication security. The baseline scheme achieves again a higher average user secrecy rate 
than the proposed scheme and simultaneous DL and UL secure transmission can be guaranteed 
when rUL < 8 dB. However, the baseline scheme consumes significantly larger DL and UL 
transmit powers than the proposed scheme, cf. Figure 151 Besides, the baseline scheme incurs 
again a high outage probability. In fact, the baseline scheme cannot guarantee communication 
security if Fj^g is larger than 8 dB since it always becomes infeasible. 

D. Average Transmit Power versus Maximum Channel Estimation Error 

In Figure lH we investigate the average transmit power consumption versus the maximum 
normalized channel estimation error, for a minimum required DL SINK of F^g = 10 dB 
and a minimum required UL SINK of Fj^g = 5 dB. The FD BS is equipped with Nt = 10 
antennas. We select again the resource allocation policy with Ai = 0.1 and A 2 = 0.9 to give 
preference to the total UL transmit power minimization. As can be observed, the average total 
DL and UL transmit powers increase with increasing maximum normalized channel estimation 
error, In fact, with increasing imperfectness of the CSI, it is more difficult for the FD BS 
to perform accurate DL beam steering. Hence, the FD BS transmits the information signal and 
the AN with higher power to satisfy the DL QoS requirements and to guarantee transmission 
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security in both DL and UL. On the other hand, the higher total DL power causes more SI 
to the UL reception. As a result, the UL users are forced to transmit with higher powers to 
satisfy the minimum required receive SINK at the FD BS. As for the baseline scheme, feasible 
solutions can only be found for < 12% which indicates that the maximum channel estimation 
error tolerance is much lower than for the proposed scheme. Besides, for the baseline scheme, 
a significantly higher power is consumed for both DL and UL transmission compared to the 
proposed scheme even if the problem is feasible. 

VI. CONCLUSIONS 

In this paper, we studied the power efficient resource allocation algorithm design for enabling 
secure MU-MIMO wireless communication with an FD BS. The algorithm design was formulated 
as a non-convex MOOP via the weighted Tchebycheff method. The proposed problem aimed at 
jointly minimizing the total DL and UL transmit powers for achieving simultaneous secure DL 
and UL transmission. The imperfectness of the CSI of the eavesdropping channels and the CCI 
channels was taken into account for robust resource allocation algorithm design. The proposed 
MOOP was solved optimally by SDP relaxation. We proved that the globally optimal solution can 
always be obtained or constructed by solving at most two convex SDP optimization problems. 
Simulation results not only revealed the trade-off between the total DL and UL transmit power 
consumption, but also confirm that the proposed FD system provides substantial power savings 
over a baseline scheme. Furthermore, the simulation results confirmed the robustness of the 
proposed scheme with respect to imperfect CSI. More importantly, our results revealed that an 
FD BS can guarantee secure UL transmission which is not possible with an HD BS. 

Appendix 


A. Proof of Proposition 1 


We start the proof by rewriting constraints C3 and C4 as 


pDL 

C4: 






det(lK,+X„'/2L«W»L„X„''2)<2 


pL 


det(I^.,+F,X-V2e^. 


,(36) 

,(37) 


respectively, (a) and {b) hold due to a basic matrix equality, namely det(I+AB) = det(I+BA). 
Then, we study a lower bound of (l3^ and (iJTl) by applying the following Lemma. 

Lemma 3 / U?]/ .- For any semidefinite matrix A ^ 0, the inequality det(I + A) > 1 + Tr(A) 
holds where equality holds if and only if Rank(A) < 1. 
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We note that 0 holds in if^ . Thus, applying Lemma 3 to t f^ yields 

det(I„, + X^'/'^LSW^L^X-'/^) > 1 + TtlX^'/'^LSW^L^X-'/^). (38) 


As a result, by combining (|3^ and (1^ . we have the following implications 

lV(X-''"L"WtL„X-‘''") < C ^ A„„(X-''"L"WtL„X->'") < 

^ X-''"L"WiL„X-‘''" ^ ^ L"WiL„ r< a°I;X„, (39) 


where Amax(A) denotes the maximum eigenvalue of matrix A and (c) is due to the fact that 
Tr(A) > Amax(A) holds for any A A 0. Besides, if Rank(Wfc) < 1, we have 

Rank(X-i/2L^WfcWX-i/2) < min I^Rank(X-^/^L^),Rank(W,L^X-^/^)j 

< Rank(WfcLyy,X-i/2) < (40^ 

Then, the equality in (1^ holds. Besides, in (1^ . Tr(Xm^^^L^WfcLmXm^^^) < are equiv¬ 
alent with Amax(Xm^'^^L^WfcLmXm^'^^) < . Therefore, (1^ and (1^ are equivalent if 

Rank(Wfc) < 1. 

As for constraint C4, we note that Rank(FjXm^'^^ej^friG^m^^^'^^) ^ ^ always holds. Therefore, 
by applying Lemma 3 to (1371) . we have 

det(I«, + P,X-‘''=e,,„e"„X-‘/=) = 1 + Tl'(PjX-‘'=ey„e"„X-‘'"). (41) 


Then, by combining (1371) and (Hdl) . we have the following implications: 


C4 Tr(F4 






^maxP 3-^m ) — Sj, 


UL 


Pe 


A ej^^Xyyy. (42) 


B. Proof of Proposition 2 


To facilitate the presentation, we define ^ and G as Sj^m = 
and Tra = ~ l-)Xm, rcspcctively. If there exists a Hermitian matrix e satisfying 

constraints C4a and C4b, then constraints C4a and C4b imply 

C4a — M 4 m)x < 0, and C4b x'^(M 4 m — Tm)x < 0, (43) 


where x G is any non-zero vector. From (l43l) . we have 


X (Sj-m - T,yj)x < 0 


C4- S ■ AT 

v_.“T. y^j^rn _ -"-m* 


(44) 
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C. Proof of Theorem 1 

The SDP relaxed version of equivalent Problem 3 in (1^ is jointly convex with respect to 
the optimization variables and satisfies the Slater’s constraint qualification. Therefore, strong 
duality holds and solving the dual problem is equivalent to solving the primal problem ll25ll . For 
obtaining the dual problem, we first need the Lagrangian function of the primal problem in 
which is given by 

K K J K 

C = Ai«iy^Tr(Wt)-y^Tr(RCT.(Wi,Z,P,,4)DcTj + 5^ft'5^Tr(pWjHgV,HE 

k=l k=\ j=l k=l 

M K K 

- EEtKR-cj,,™ ki Z, ^ )-5^Tr(WfeYfc) + A. (45) 

m=l k=l k=l 

Here, A denotes the collection of terms that only involve variables that are independent of W^. Hj 
and Oi are the Lagrange multipliers associated with constraints C2 and C9, respectively. Matrices 
Dot, 6 €<"+')=< (>^+ 1 ) and G C^^t+^r)x(^t+^^r) are the Lagrange multiplier matrices for 

constraints Cl and C3, respectively. Matrix G is the Lagrange multiplier matrix for 

the positive semidefinite constraint C7 on matrix W^. For notational simplicity, we define as 
the set of scalar Lagrange multipliers for constraints C5, C9, and CIO and as the set of matrix 
Lagrange multipliers for constraints Cl, C2, C3, C4a, C4b, C 6 , and C7. Thus, the dual problem 
for the SDP relaxed problem in (l35l) is given by 


(46) 




maximize minimize £(Wfc, Z, M^-r, 4,4^^, T', # 


s.t. 


k, ■ 

to- 

2 = 1 


= 1 . 


Constraint Ytd=i 6*1 = 1 is imposed to guarantee a bounded solution of the dual problem 
Then, we reveal the structure of the optimal of (l35l) by studying the Karush-Kuhn-Tucker 
(KKT) conditions. The KKT conditions for the optimal are given by: 

(47) 

Y*W* = 0, (48) 

Vw‘/: = 0, (49) 

where Yl, D^r , , /U*, and 6* are the optimal Lagrange multipliers for dual problem (l46l) . 

LaI/c kyOk,m ^ 

denotes the gradient of Lagrangian function C with respect to matrix W^. The KKT 
condition in (l49l) can be expressed as 


M 


^ diag(V,)Hsi + Bl„ 


D4, 


D2 


C3k,m riH I -T) Glk rtH 


i=i 


m=l 


pDL 

reqfc 


pDL hfc- 
reqfc 


(50) 
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Now, we divide the proof into two oases acoording to the value of Ai. First, for the ease of 


0 < Ai < 1, we define 
,/ 


M 


D— 

C3fc, 


AI = j;,.>HSdiag(Vj)Hs, + and 


J = 1 


m=l 


for notational simplioity. Then, (1501) implies 

Y* = nt - 


D— 

^hfc ynDL 


(51) 


(52) 


Premultiplying both sides of (1521) by W^, and utilizing (14^ . we have 


W*n* = 


k -pDL hfc • 


(53) 


By applying basio inequalities for the rank of matrioes, the following relation holds: 


Rank(W*) = Rank(W*n*) = Rank (W*B 


D— 

Clk -oH 
-pj^L 
reqfc 


(h) ( / \ 'i (c) / D7 

< mm|Raiik(Wa.Rank(BH. 7 ;i|iB»)} < Rank .(54) 


^B«)}<Rank(Ba,^-" 

reqfc ^ req*. 

where (a) is due to 11^ 0, (&) is due to the basio result Rank( AB) < min | Rank( A), Rank(B)}, 

and (c) is due to the faot that min{a, b} < a. In order to further reveal the struoture of in 

D:^ 

, we study the rank of Bh^. pn^B^. whioh is given by 
Rank ( Bh, .^^B^ ) = Rank ( [O h,] 


pDL 

reqj. 


Clfc 'qH 


pDL hfe 

(d) D— ’'''le) 

< min I Rank h^] j, Rank | < Rank [O h^] j < 1, (55) 

reqfc 

where for (d) and (e), we used the same results as for (6) and (c), respeotively. By oombining 

(56) 


(l54l) and (l55l) . the rank of is given by pj., 

Rank(W:) < Rank(Ba,p^B") <1. 


reqfe 


We note that ^ 0 for F^^^ > 0. Thus, Rank(W^) = 1. 

Then, for the ease of Ai = 0, we show that we oan always oonstruet a rank-one optimal 
solution W". We note that the problem in (1351) with Ai = 0 is equivalent to a total UL transmit 


power minimization problem whioh is given by 

minimize \ P, 

Wfc,ZeH^T,Mi 

J. 


s.t. 


O, C2, C3, C4, C5, C6, C7, C9, CIO. 


(57) 


We first solve the above oonvex optimization problem and obtain the UL transmit power P**, the 
DL beamformimg matrix W^, the AN eovarianee matrix Z*, and the optimal auxiliary variables 
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which are collected in H* = {M*^, 61,1% ,^, If Rank(W^) = 1, Vfc, then the globally 

optimal solution of problem (|2T]) for Ai = 0 is achieved. Otherwise, we substitute P**, Z*, and 
S* into the following auxiliary problem: 

minimize y^Tr(Wfc) + Tr(Z*) 

WfcGH^T ^ 

s.t. O,C2,C3,W,C5,C6,C7,C9,C10. (58) 


Since the problem in (1581) shares the feasible set of problem (1571) . problem (l58l) is also feasible. 


Now, we claim that for a given P**, Z*, and H* in (158]) . the solution W” of (1581) is a rank- 


one matrix. First, the gradient of the Lagrangian function for 
expressed as 




Bh 


Clk-p,H 

Bh;,) 


k -pDL 

reqfe 


with respect to W^* can be 

(59) 


where n;* = + AJ* and AJ* = E"=i diaglWjHs, + E".i Bl 


and 


Yr, , and /j,** are the optimal Lagrange multipliers for the dual problem of 


Premultiplying both sides of 


by the optimal solution W^*, we have 

jy** 

Clfc pg 

pDL ^h;,' 


wrnr = wrBh, 


(60) 


reqfe 


We note that 11” is a full-rank matrix, i.e., 11” 0, and (1^ has the same form as (l53l) . Thus, 

we can follow the same approach as for the case of 0 < A* < 1 for showing that W” is a 
rank-one matrix. Also, since W” is a feasible solution of (1571) for P”, an optimal rank-one 
matrix W” for the case of Ai = 0 is constructed. 
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